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SUMMARY 

The lethargy-dependent equations of the consistent PI approximation 
to the Boltzmann transport equation for slowing down neutrons have been 
used as the basis of an IBM 704 computer program. Some of the effects 
included are (l) linearly anisotropic center of mass elastic scattering, 
(2) heavy element inelastic scattering based on the evaporation model of 
the nucleus, and (3) optional variation of the buckling with lethargy. 

The microscopic cross-section data developed for this program covered 
473 lethargy points from lethargy u - 0 (10 Mev) to u *= 19.8 (0.025 ev). 

The value of the fission neutron age in water calculated here is 
26.5 square centimeters; this value is to be compared with the recent 
experimental value given as 27.86 square centimeters. The Fourier trans- 
form of the slowing -down kernel for water to indium resonance energy 
calculated here compared well with the Fourier transform of the kernel 
for water as measured by Hill, Roberts, and Pitch. 

This method of calculation has been applied to uranyl fluoride — 
water solution critical assemblies. Theoretical results established for 
both unreflected and fully reflected critical assemblies have been com- 
pared with available experimental data. The theoretical buckling curve 
derived as a function of the hydrogen to uranium— 235 atom concentration 
for an energy— independent extrapolation distance was successful in pre- 
dicting the critical heights of various unreflected cylindrical assem- 
blies. The critical dimensions of fully water-reflected cylindrical 
assemblies were reasonably well predicted using the theoretical buckling 
curve and reflector savings for equivalent spherical assemblies. 


INTRODUCTION 

The NASA Lewis Research Center has for some t im e operated a Zerp 
Power Reactor. This reactor consists of uranyl fluoride salt (U0 2 F 2 ) 
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dissolved in water and contained in a cylindrical aluminum tank. The 
parameters of the system are the diameter of the cylindrical tank, use 
of a side reflector, and concentration of uranyl fluoride salt in water. 
Criticality is achieved by varying the height of the solution in the 
tank. No control rods are inserted into the solution while the reactor 
is critical. Thus, good experimental conditions exist for obtaining 
criticality data and performing reactivity measurements by measuring 
changes in the solution height due to insertions in the reactor. Descrip- 
tion of this reactor and its various characteristics are found in the 
final hazards report (ref. l) prepared for the Atomic Energy Commission. 

The criticality data obtained from the NASA Zero Power Reactor as 
well as that obtained from extensive experimental studies performed at 
the Oak Ridge National. Laboratory (refs. 2, 3, and 4) can be analyzed 
by various theoretical models. Callihan (ref. 5) describes the use of 
Fermi age-diffusion and Goertzel-Selengut analytical models as well as 
a method due to Greuling for predicting critical reactor diameters of 
various Oak Ridge experiments. A more recent study of the Oak Ridge 
experimental data by Gwin and others (ref. 6) makes use of the empirical 
kernel method, which is based on the measured slowing-down distribution 
in water rather than on any theoretical model of the neutron slowing- 
down process. In addition, the use of the Los Alamos DSN transport code 
for obtaining minimum critical dimensions for water solution is described 
by Mills (ref. 7). 

The objective of the present study is to provide a theoretical model 
based on the numerical solution of the lethargy-dependent consistent PI 
approximation to the Boltzmann transport equation for slowing down 
neutrons. This model is coded for the IBM 704 electronic computer and 
employs microscopic nuclear data (refs. 8, 9, and 10). This theoretical 
model is tested by comparison with the experimental age and Fourier 
transform of the slowing-down kernel of fission neutrons in water. This 
model is further employed in estimating the absolute reactivity of un- 
reflected UO2F2 -wat er systems and is used as a basis for evaluating con- 
sistent three-group parameters for the interpretation of water-reflected 
UOfFf -water systems. 

The next section gives the derivation of the general lethargy- 
dependent consistent PI equations used as the basis for an IBM 704 code. 
Next follows a section giving the analytical results based on this model 
and comparison with experimental data obtained with the NASA Zero Power 
Reactor as well as data obtained at Oak Ridge. 


THEORY 

Slab Geometry Spherical Harmonics Form of Transport Equation 
> — *■» - > —*■ 

Let n(r, vjdr dv be the probable number of neutrons in the space 
volume element dr about position r and in the velocity element dv 
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about velocity v per unit time. Then, the Boltzmann transport equa- 
tion for slowing down neutrons at steady-state, which equates the losses 
per unit time of neutrons, due to transport out of and collision re- 
moval from the phase volume element d? dv, to the gain per unit time 
of neutrons, due to elastic and inelastic in-scattering and fissioning 
of the fissile nuclei in the phase volume element dr dv, has the form 

n • Vq)(?,^) + = f dv'£ es (v' )cp(r, 7’ )P(?< -» ?) 


+/i?'S ln (T> )<p(?, V<)H(?' -?) + S(?,v) (1) 

where the neutron flux is defined as cp(?,v) = vn(r,v) with v being 
the neutron speed; S(r, v) is the source of neutrons due to fission; 

^t( v )j £ es (v), £ in (v) are the macroscopic total, elastic scattering, 

and inelastic scattering cross sections, respectively, at speed v; Q 
is a unit vector along the direction of motion of the neutron; Pfv* v) 

is the probability that a neutron in the phase volume element dr dv r 
w i i:L 4 after an elastic scattering collision, be in the phase volume 
dr dv. Similarly, H(v’ -*> v) is the probability that a neutron in the 
phase volume element dr dv f will, after an inelastic scattering col- 
lision, be in the phase volume dr dv. 

Now t ran s f orm^equat i on (l) from (r,v) space to (r,v,n) space and 
finally from (r, v, Cl) space to (r, u, {}) space where u is the neutron 
lethargy defined by the relation u = ln(E Q /E) with Eq usually taken 

as 10^ ev. The result of these transformations and writing an explicit 
form for the fission source is 

ft * Vcp(r, u, ft) + £t(u)cp(r,u,ft) =f d.u'2 es (u') f dfi’ cp(?, u;ft' )P(u-u; n L ) 

+ /du«s in (u>) H ( u ' u ) f cm’tp(r,u;fl> ) 

+ lS:/ du, ^ u ’ )s fC u ')/m'9(7,u;n') ( 2 ) 

where v(u) is the number of neutrons released per fission occurring at 
lethargy u, Zf(u) is the macroscopic fission cross section, f(u) is the 

fraction of the fission neutrons born in the lethargy element du about 
u (the integral of f(u) over all lethargies is normalized to unity), 

Kgff is the static criticality factor, and dn is the element of solid 

angle about the direction ft. The inelastic scattering kernel is assumed 
to depend only on the initial lethargy u* and the final lethargy u 
of an inelastically scattered neutron and to be isotropic in the labora- 
tory coordinate system. The elastic scattering kernel is assumed to 


4 


depend only on the difference bet-ween the final and initial lethargies 
and on the cosine of the angle jj.^ between the initial and final di- 
rections of motion of a neutron undergoing an elastic scattering col- 
lision. Davison (ch. IX, ref. ll) gives additional details concerning 
equat ion (2). 

Xf the neutron distribution varies only along the Z-direction and 
is azimuthally symmetric about that axis,, if equation (2) is integrated 
over the azimuthal angle, and if the flux cp (- cp(z, u, q) (the quantity 
(j. is the cosine of_^the angle between the Z-direction and the direction 
of neutron motion ft) and the elastic scattering kernel are expanded 
into a series of Legendre polynomials, that is, 


uu 

cp(Z,u,u) = £ C p J (z,u)Pj(|a) 

2=0 

(3a) 

00 

P(u-U,V L ) = E 2 1 ^(u-u')Pj(n L ) 

(3b) 


2= 0 


and these expansions are substituted into the equation for cp(z, u, q), 
which is then multiplied by a Legendre polynomial P K (q) and integrated 

over all q, then the result will be the infinite sequence of integro- 
differential equations given as 

^zfjfT+T + 2 y + 1 ^2+l( Z ' u )] + ^feM^C^u) 

=/ du,£ es(u f )cpz(Z,u* )P z (u-u t ) + S 0 ,1 f duf ^in( uf )9o(^ )h(u t - u) 

+ dU ' V(U ' )£f(U ' )cp 0^ u ’)^ (1 ^ 0, 1, Z, (4) 

where 5 Q j is the Kronecker delta defined here by the properties 

5 0, Z = 1 = 0 

= 0 , 2 £ 0 ( 5 ) 

Equation (4) represents the spherical harmonics form of the transport 
equation in slab geometry. 
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Slab Geometry Spherical Harmonics Form of Angular Slowing- 
Down Density Due to Elastic Scattering 

Another important characteristic of slowing down neutrons is the 
slowing-down density Q(r,u,jf). The slowing -down density due to elastic 
scattering Q(r, u, rf) is defined as the number of neutrons which pass the 
lethargy mark u per unit time per unit phase volume (dr dfi). This 
quantity is easily derivable as 




/ dG 'Z es (u f )cp(r, uj )P(u* 


( 6 ) 


If the slowing-down density varies only along the Z-direction and 
is azlmuthally symmetric about that axis, if equation (6) is integrated 
over the azimuth angle cp, if the expressions for cp(Z, u, q) and 
P(u-u’n L ), given by equations (3a) and (3b), respectively, are sub- 
stituted into (6) and then integrated over 1 , and if the resulting 
equation is multiplied by a Legendre polynomial P K (q) and integrated 

over all q, then the result will be the moments (Z, u) of the 
slowing-down density Q(Z, u,?f); that is, 

^U oo 

Qj(Z,u) = J dnQ(Z,u,jf)Pj(n) = J du J du"£ es (u’ )<p, (Z, u' )Fj (u" - u’ ) 

^ 0 *ui 


(2 = 0 , 1 , 2, . . . , oo ) ( 7 ) 

Equation (7) represents the spherical harmonics form of the slowing-down 
density in slab geometry. 


Lethargy- Dependent Consistent PI Equations in Slab Geometry 

To obtain the PI equations from equation (4) the assumption is now 
made that cp^(z, u) = 0 for Z ^ 2, which implies for equation (7) that 
^2 (Z-> u) = 0 for 2^2 also. The lethargy-dependent consistent PI 
equations can now be obtained by assuming the following spatial and 


lethargy dependence for cp Q (Z, u) and cp 1 (Z, u): 

cp Q (Z,u) = <p(u)G[B m (u)z] 

(8a) 

cp-L^u) = j(u)F[B m (u)z] 

(Sb). 
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In addition, the material buckling B^u) is assumed to be a slowly 
varying function of lethargy. Therefore, the integral terms arising in 
equation (4) can be evaluated by taking 


<Po( z > u ') = 

<Pl(Z,u') & j(u')?[b,Wz] 


(8c) 

(3d) 


The lethargy and space variation for the l - 0,1 terns of equation (7) 
follows from equations (8a) and (8b) and is 


Qq(Z,u) = Q(u)G[B m (u)z] 
Q 1 (Z,u) = P(u)F[B m (u)Z] 


(8e) 

(8f) 


where G and F are spatial functions and B^u) is defined as the 
lethargy -dependent material buckling. In addition, the spatial functions 
G and F must satisfy the conditions 


52 {o[Bm(“)Z]} = " BmWFOBmCuJz] 
i {fK(u)z]} = B m (u)G[E m (u)z] 


(9a) 


(9*>) 


With the assumptions implied by equations (8a) to (8d) and (9a) and (9b), 
the spatial dependence can be divided out, and the lethargy-dependent 
consistent PI equations obtained from equation (4) are 

B m (u)j(u) + Z t (u)cp(u) = / clu'S^u* )<p(u' )P 0 (u-u’) 


+ J m du r E in (u ? )cp(u f )h(u' -► u) 

+ f(u| r (i u » v ^ u «j2 f(u')cp(u T ) 
■^eff Jq 


BjufaMu) + Z t (u)j(u) =J du , Z es (u f )J(u , )i > l(^“■w , ) 


(10a) 

(10b) 


together with 


U co 

Q(u) = J du f J du ,, Z es (u' )cp( uf )Po( u " ~ ) (Ha) 


u 


069-a 
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/ U oo 

•3-V.' j du"Z es (u')J(u')Pi(u" - u') 


(lib) 


In reference If Simon discusses some analytical solutions of the 
lethargy-dependent consistent PI equations for hydrogen. To obtain the 
analytical results, special cross-section variations were assumed for 
hydrogen. 


Elastic Scattering Kernel for Linearly Anisotropic 
Center of Mass Elastic Scattering 

Since most light elements exhibit anisotropic elastic scattering 
in the center of mass coordinate system for high-energy neutrons, it is 
important to include this effect if accurate values of the neutron age 
are to be computed. Therefore, by following Marshak (ref. 13) it is 
easy to show that for linearly anisotropic scattering in the center of 
mass coordinate system the elastic scattering kernel P(u-u'p T ) is 
given by L 

P(u-u,Vp)du dp^ dcp-^ 



X du dp L dcf^ ( 12 ) 


where fq(u) is the average cosine of the angle between the initial and 
final direction of motion in the center of mass coordinate system of the 
elastically scattered neutron of initial lethargy u, A is the ratio 
of the mass of the scattering nucleus to that of the neutron, 5 is the 
Dirac delta function, and is an appropriate azimuthal angle in the 

laboratory coordinate system. The moments p 2 (u— u f ) of the expansion 
of P(u-u, Pp) into a series of Legendre polynomials, as in equation (3b), 
are given by the orthogonality property of these polynomials as 

P z (u-u') =y*dn L P z (p L )p(u-u, f p L ) (13) 

Substituting equation (12) into equation (13) and performing the indi- 
cated integrations for l = 0, 1 will give 
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Pq(u-u') = ( A e - (u - U ’ 5 |l + 3f- L (u') [( A l A l)2 e -( u - u ') 


A 2 + 1 
2A 


>, -1 5 M-t (u-u* ) < 1 


- 0, otherwise 

P-^u-u’ ) = )P q (u-u t ), -1 < p. L (u-u f ) 5 1 

= 0, otherwise 


(14a) 


(14b) 


If t = u-u f , equations (l4a) and (l4b) can he rewritten into the 
following useful forms: 


P n (t) = e- 6 + e -2t - 3(A 2 4l)(A + l) 2 e _ t l 

u ^ L 8A 2 8A 2 J 


(15a) 


3 

ft) - ( A + e “2 * . (A - 1)(A + 1)* “2 

\ Z} ~ 8A e «a e 


2 4t 


8A 


+ f-^u* ) 


•5(A + l) 5 “f 13 _ 5(A + l) 5 “| t 


16A^ 


8 


1 

2 4t 


+ 


5 (A - l) (A 2 + l) (A + l) 2 ~2 


16A Z 


(15b) 


Equations (15a) and (15b) reduce to their isotropic center of mass 
elastic scattering values if f i( u ' ) s °- 


Differential Form of Consistent PI Equations 

Although equations (10a) and (10b) are the lethargy-dependent con- 
sistent PI equations, a differential form of these equations would be 
more suitable for computation. !fo obtain these differential equations 
differentiate equations (lla) and (lib) with respect to lethargy u: 

~ y U ) p oO i ** u, ) + S es ( U M U ) ./* du ,f P 0 (u" - u) 

U.U 0 U 

= - S* U ^u , Z es (u' )cp(u T )P 0 (u-u T ) + s (u)cp(u) (16a) 

0 
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- f 0 du'2 es ( U ')J(ti')P 1 (u-u') + S es ( u )j(u )jf du" Pl (u" - u) 

= "y au ' z es( u ') J (' u ') p 1 ( u - l;l ') + |a 0 (u)£ es (u)j(u) (16b) 


The integral J du"P 0 (u" - u) is evaluated using equation (15a) and 

U ^ y^OO V 

is equal to unity. The integral / du ,, P 1 (u M - u) is evaluated using 

*ui 

(l5h) and is equal to |a 0 ( u ) “ average cosine of the angle between 

the initial and final direction of motion, as measured in the laboratory 
coordinate system, of the elastically scattered neutrons. The first 
integral on the right side of equation (lOa) can be evaluated using equa- 
tion (l6a), and the integral term on the right side of equation (10b) 
can be evaluated using equation (l6b). Equations (lOa) and (lOb) now 
become 


B ra (u)j(u) + [E a (u) + Z ln (u)]cp(u) + MHl 

= fid f du'v(u')S f (u')<p(u') + /’du'S ln (u')cp(u')H(u' - u) 

Iv eff *y) 

(17a) 

- 3 Bm( u )cp(u) + Pt( u ) " P-o( u ) E es( u )] J (") + ^d = 0 (17b) 


where £ a (u) is the macroscopic absorption cross section. Equations 

(17a) and (17b) represent the basic lethargy-dependent equations of the 
consistent PI approximation to the Boltzmann transport equation for 
slowing -down neutrons. 


Two more equations are needed to complete the PI description of 
the slowing-down process as given by equations (l7a) and (lTb). These 
are coupling equations between cp(u) and Q(u), and J(u) and P(u), 
which can be obtained by writing approximate forms for equations (ll) 
and (16). These additional equations are obtained in the following 
manner. Assume that the functions Z es (u 1 )<p(u' ) and E es (u r )«T(u r ) are 

slowly varying functions and that they can be approximated by the first 
two terms of a Taylor's series expansion about lethargy u. That is, 

S es (u') 9 (u') i E es (u)<p(u) + (u' - u) A [S es (u)cp(u)] (18a) 

SesCuOjfu 1 ) = S es (u)j(u) + (u 1 - u) d. p es (u)j(u)] (18b) 
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Substitute equations (l8a) into (lla) and (l6a) and equation (l8b) into 
(lit)) and (l6"b). Then,, evaluating the various integrals over lethargy 
of Pq(u-u’ ) and P-^u-u' ) and eliminating the terms involving 

C^es ( U ) C P( U )] and — - [S es (u)j(u)] result in 

r (u) + Q(u) = l 1 (u)S es (u)cp(u) (19a) 

p(u) + P(u) = q 1 (u)Z es (u)j(u) (l9h) 

The parameters appearing in equations (19a) and (19b) are defined, in 
terms of the quantity t = u-u T , as averages over the distribution 
function Po(t). That is, 


^1 - t 

(20a) 

2£-jj - t2 

(20b) 

q s tq L (t) 

(20c) 

qP = t 2 ^ L (t) 

(20d) 


where the bar denotes an average value. Of course, the various param- 
eters defined by equation (20 ) depend on lethargy since Pq(u-u t ) and 
P^u-u’ ) depend upon fq(u). 

Thus, equations (l7a), (l7b), (l9a), and (l9b) form the four 

coupled first-order ordinary differential equations describing the 
neutron slowing-down process in the consistent PI approximation to the 
Boltzmann transport equation. Included are the effects of linearly 
anisotropic center of mass elastic scattering, inelastic neutron scat- 
tering, and lethargy-dependent geometric buckling. Although equations 
(17) and (19) apply for only one nuclear specie, the extension of the 
equations to include additional nuclear species follows directly. In 
particular a set of equations (19) is written for each element with 
appropriate modifications of equations (17). 


Comparison of Lethargy-Dependent Consistent PI Equations 
With Various Approximate Theories 

Equations (17) and (19) are fairly general, and it is worth noting 
that under various restrictions simpler sets of equations, which have 
been used frequently in reactor analysis, can be easily obtained. This 
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will "be done for a few typical cases. The first of these is the age- 
diffusion equation with no inelastic scattering which can he obtained 
from (17) and (19) by setting £ in (u) =0 in (l7a), dP(u)/du =0 in 

(17b), dQ(u)/du =0 in (l9a), and dP(u)/du = P(u) =0 in equation 
(19b). 

The familiar Fermi age-diffusion equations result and are given as 
(ref. 14) 

B m (u)j(u) + E a (u)cp(u) + £ du'v(u')2 f (u')cp(u') (21a) 

J(u) = D(u)B m (u)cp(u) (21b) 

Q(u) = i 1 (u)E es (u)cp(u) (21c) 

where the diffusion coefficient D(u) has been defined as 

^ ^ 3 ta( u ) + S es( u )C 1 ~ M- 0 ( u )]} < ’ 2Z) 

The Selengut-Goertzel approximation is basically the age-diffusion set 
given by equations (21a) to (21c) but with the element hydrogen treated 
by equation (19a) - other elements still being treated by equation (21c). 
To see this, recall that hydrogen exhibits isotropic center of mass 
elastic scattering collisions over energies of interest in reactors. 

Thus, £p(u) and y(u) for hydrogen are both equal to unity for all 
lethargies giving 

+ Q H ( u ) = sf s (u)cp(u) (23) 

which is the Selengut-Goertzel equation relating the hydrogen slowing- 
down density Q H (u) to the lethargy-dependent flux cp(u) (p. 348, ref. 
15). The Greuling-Goertzel approximation (p. 390, ref. 15) includes 
both equations (21a) and (21b) of the age-diffusion formalism but re- 
tains equation (19a) for all elements. Other approximations are also 
possible but will not be discussed because machine programs are just as 
easily formulated for the basic set of equations given by (17) and (19 ) 
as for even the simple age -diffusion set of equations. 


Equations for Thermal Group 

The slowing-down equations given by (l7a), (l7b), (19a), and (19b) 
are assumed to apply to some cutoff lethargy up corresponding to an 

energy Ep. The energy Ep is usually taken as about 5kT where T 

is the temperature of the moderator in °K and k is the Boltzmann con- 
stant. In the energy range from 0 to Ep the neutrons have energies 
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comparable to the energies of the thermally agitated nuclei and can 
gain as well as lose energy upon making a collision. However, if the 
absorption and leakage of neutrons are small, the neutron energy spec- 
trum will approach a Maxwell-Bolt zmann distribution in speed with a 
temperature corresponding to that of the moderator. Under these con- 
ditions and using the principle of detailed balance, the thermal equa- 
tion in the PI approximation will be, after separating the spatial 
function as before, 

du / du r 2 in (u’ )cp(u* )H(u T -» u) (24a) 

*^-00 

*3 + ^trrp'~*T “ ^’( u r) (24b) 

where and are the appropriate averaged values of the macro- 

scopic absorption and transport cross sections over the thermal speed 
range. Thus, equations (17a), (17b), (19a), (19b), (24a), and (24b) 
form the complete set of consistent PI equations extended to include 
many elements, which are the basis for an IBM 704 code. 


B m T J T + 2 a^ ~ 


+ 


jL 


b= 


* 


Discussion of Boundary Conditions, Fission Spectrum, 
and Inelastic Scattering Kernel 


If the zero of lethargy is taken sufficiently lower than the low- 
est lethargy in the source spectrum, the initial conditions for the 
set of consistent PI equations will be 

Q(u=0) = 0, cp(u=0) = 0 (25a) 

P(u=0) = 0, J(u=0) = 0 . (25b) 


Depending upon the reactor, a number of possible source spectrums can 
be used. For this study the uranium-235 fission spectrum, given in 
analytic form as 


f (u) 


= 0.45270 E 



(26) 


has been used. The quantity E is the neutron energy in Mev. The 
spectrum is normalized to unity; that is, 




(27) 


This normalization of the source spectrum permits the evaluation of the 
static criticality factor of the reactor as 


;-byu 
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/ CO 

duv(u)S f (u)cp(u) (28) 

The characteristic dimension of the reactor is needed to compute the 
lethargy-dependent buckling B^u). For example,, in the case of a 
spherical reactor the equation used for B^u) is 


TL(u) = r— - 

m r + 0.71 yu) 

■where r is the radius and "where 

\ r (u) = 3D(u) 


(29) 


(30) 


The inelastic scattering kernel in equation (l7a) has "been obtained 
from the evaporation model of the nucleus as described by Blatt and 
Weisskopf (ref. 16). The evaporation model of the nucleus appears to 
give the gross features of heavy element inelastic scattering. In terms 
of energy, this model gives for the inelastic kernel the following 
result: 

H(E r - E) = const. E[e"^ E ^ E ] (3l) 

where the constant is determined from the normalization requirement that 


J dE H(E* - E) =1.0 (32) 

In equation (31) E r represents the initial energy of the neutron and 
E the energy after an inelastic scattering collision has occurred. The 

quantity £(E’) is the reciprocal nuclear temperature in Mev~-^- and is a 
function of the scattering nucleus and initial energy of the neutron. 

It has been approximated by 

£(*')- (i§r) l/2 (33) 

with the parameter a being experimentally determined for each nucleus 
of interest. 


Discussion of Neutron Cross Sections 

The lethargy-dependent consistent PI equations described along 
with the initial conditions, source spectrum, and approximate inelastic 
scattering kernel are integrated for a homogeneous unreflected reactor 
system over the microscopic cross-section data of the elements involved 
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in a stepwise fashion from 10 7 ev to 0.025 ev (or from lethargy u = 0 
to u - 19.8). This lethargy range is spanned by 473 mesh points such 
that the cross sections of the various elements are well represented 
in the resonance energy range. The microscopic cross-section data of 
29 elements are represented in such a fashion for use with the IBM 704 
code. These cross-section data include a a , a gg , vo f ljr and o ±n 

where cr a is the microscopic absorption cross section, a es is the 

microscopic elastic scattering cross section, va^ is the microscopic 

production cross section, f-j_ is the average cosine of the angle between 

the initial and final direction of motion of an elastically scattered 
neutron as measured in the center of mass coordinate system, and ov 

is the microscopic inelastic scattering cross section. The quantity v 
is the average number of neutrons produced in fission, and o ^ is the 

microscopic fission cross section. The quantity v is taken to be 
energy-dependent according to the data of reference 17. In the thermal 
energy range the absorption and fission cross sections are corrected for 
non-l/v behavior as well as for hardening of the thermal spectrum, due 
to absorption, using the data of reference 18. The transport cross sec- 
tion at thermal energy is similarly corrected. 


THEORETICAL RESULTS AND COMPARISON WITH EXPERIMENT 
Age of Fission Neutrons in Water 

One of the important parameters describing moderating materials is 
the fission neutron age t. For hydrogenous media the age becomes espe- 
cially important because of the well-known discrepancy for water between 
the experimental measurements and various theoretical calculations. For 
any homogeneous medium the lethargy-dependent consistent PI equations 
can be solved for a given set of microscopic cross-section data, the 
uranium-235 fission spectrum, and for a constant buckling Bj^. From 
the neutron fluxes and currents generated as a function of lethargy, it 
is possible to compute the weighted fission age. 

In this work, the weighted fission age for water to indium reso- 
nance energy is computed to be 26.5 square centimeters. Oxygen was con- 
sidered to have linearly anisotropic center of mass elastic scattering. 
For isotropic center of mass scattering for oxygen, the age is calcu- 
lated to be 24. 7 square centimeters. This demonstrates the importance 
of precise calculation of the various competing effects on the neutron 
spectrum. For comparison of the age of 26.5 square centimeters calcu- 
lated in this work, the weighted mean of a number of different types of 
calculations is given in reference 19 as 26.Q±0.5 square centimeters. 

It is also of interest to note that a recent measurement of the fission 
neutron age in water is 27.86 square centimeters (ref. 20 ). 


E-690 



E-S90 


15 


Fourier Transform of Slowing -Down Kernel in Water 

Another important property of homogeneous^ moderating media can he 
described in terms of the Fourier transform K^) of the slowing-down 

kernel. Physically, the Fourier transform K(B m ) is the nonleakage 

probability of the source neutrons to a given energy. The transform 
■^(®m) can derived from a solution of the lethargy- dependent consistent 
PI equations for a given source spectrum and a constant value of the 
buckling B^. 

Figure 1 is a plot of the transform K(B m ) as a function of ^ 
for fission neutrons in water to indium resonance energy. Shown on the 
curve are values of the transform computed by Gwin, Trubey, and Weinberg 
(fig. 1, ref. 6) from the experimental slowing-down kernel as measured 
by Hill, ^ Roberts, and Fitch (ref. 2l). The agreement over a wide range 
of buckling is seen to be good. The value of the fission neutron age in 
water reported in reference 21 is higher than the recent value of ref- 
erence 20 by about 3 square centimeters. This is believed due to dif- 
ferences between the two experiments of measurements near the source. 

This difference near the source of the reference 21 measurement of the 
slowing-down kernel in water is expected to have a small effect on the 
calculation of the transform ^(Bq). 

Table X compares the transform K(B rn ) to indium resonance energy 
for water and two concentrations of uranyl fluoride salt in water solu- 
tion. The differences as compared to water are slight except for extreme 
values of the buckling. 


Calculation of Critical Dimensions of Unreflected 
Homogeneous UC^Fg-HgO Assemblies 

The lethargy-dependent equations of the consistent PI slowing-down 
model described in the section on theory will be used to calculate the 
parameters of critical assemblies of unreflected uranyl fluoride salt 
UO 2 F 2 dissolved in water. Comparison of theoretical results with avail- 
able experimental data will also be made. 

Before any calculations can be performed, the physical properties 
of uranyl fluoride — water solutions must be known. Figure 2 is a curve 
of density of uranyl fluoride - water solutions as a function of hydrogen 
atom concentration for a temperature of 25° C obtained from available 
experimental data (ref. 22). The point on the curve for pure water is 
for a density of 0.997 gram per cubic centimeter and a hydrogen atom 
concentration of 6.66X10 atoms per cubic centimeter. Since the uranium 
in the solution is considered to contain 93.2 percent by weight of the 



uranium-235 isotope, the data of figure 2 can be used to derive a curve, 
figure 3, of the uranium-235 atom concentrations R 28 in atoms per 
cubic centimeter as a function of the hydrogen atom concentration R^ 
in atoms per cubic centimeter. On figure 3 are also indicated various 
values of R defined to be the ratio of the hydrogen atom concentration 
to uranium-235 atom concentration, R — R-^/r 28 . Various important char- 
acteristics of the UO 2 F 2 -H 2 O critical assemblies will be given in terms 
of this parameter R. 

Thus, specification of the solution density for the U02^2~-^2 8 solu- 
tion reactor in question gives the hydrogen atom concentration from the 
curve given by figure 2. Then, from figure 3 the uranium-235 concen- 
tration can be determined. The atom densities for the other constituents 
of the solution are determined from the following relations: 


Uranium- 2 38: 

: R 28 - 0.0730 R 28 

(34a) 

Fluorine: 

RF = 2(R 25 + R 28 ) 

(34b) 

Oxygen: 

R° 2 - R H + K F 

(34c) 


By using the atom concentration for the various constitutents of 
the UO 2 F 2 -H 2 O solutions in the manner described, calculations were per- 
formed, using the lethargy-dependent consistent PI equations over the 
473 lethargy points of the microscopic cross-section data extending from 

10 ev to 0.025 ev, over a wide range of solution concentrations. In 
terms of the parameter R, the solution concentrations were varied from 
R = 25 to 700. As a result of these many calculations, the static 
criticality factor K ef f is obtained for a given hydrogen to uranium- 
235 atom concentration ”(R) as a function of the buckling B^ of the 
system. The buckling, although capable of being varied with the lethargy 
in the IBM 704 code, was considered to be lethargy- independent for these 
calculations. 

As examples, the lethargy-dependent flux cp(u) is plotted as a 
function of lethargy u for critical reactors having solution con- 
centration values of R = 40 and 472. The source term for the cal- 
culations corresponds to 1 fission neutron. The flux for R = 40, a 
system heavily loaded with uranium-235, is shown in figure 4. In the 
fission spectrum lethargy range from about lethargy u - 0 to u = 6 
a resonance structure in the flux is noted which corresponds to various 
high-energy scattering resonances of oxygen. In the resonance energy 
range from lethargy u = 6 to u = 18.2, the decrease in the flux with 
lethargy occurs because of the absorption by the uranium-235 nuclei. 

The dips in the flux correspond to absorption resonances in the uranium- 
235. The effect on the flux of the approximately l/v capture is 
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evident from lethargy u = 15.0 to 18.2. The flux for R = 472, a 
system lightly loaded with uranium-235, is shown in figure 5. In the 
lethargy range from u — 0 to 6 the flux for R = 472 appears roughly 
the same as for the R = 40 system. In the lethargy range from u — 6 
to u = 18.2 the flux is constant from about lethargy u = 7 to 
u = 12.0. Pram lethargy u = 12.0 to 18.2 the effects of uranium-235 
absorption resonances and l/v capture are evident but are not as 
drastic as for the heavily loaded R — 40 reactor system. 

Figure 6 is a plot of the static criticality factor K eff for un- 
reflected U02F2-H20 assemblies as a function of the parametric R and 
the buckling B^. Curves of the criticality factor varying from 0.98 
to 1.12 are shown. All the curves exhibit a broad minimum in the buck- 
ling over values of R from 40 to 70. For the static criticality 
factor equal to unity the radius of the equivalent unreflected sphere, 
which includes the extrapolation distance, varies from about 18. 25 to 
25 centimeters over the range of the parameter R. 


Comparison of Calculated Critical Dimensions with Experimental 
Values for Unreflected U0 2 F 2 -H 2 0 Assemblies 

In order to compare calculated critical dimensions with experi- 
mentally determined values, the effect of the aluminum containment 
vessels on measured critical dimensions of unreflected U0 2 F 2 -H 2 0 assem- 
blies must be known. The worth of the aluminum bottom plate for the 
NASA experiments was measured in terms of the change in critical height 
as a function of the thickness of the bottom plate for a 12-inch- 
diameter cylindrical tank. The solution concentration R was 475. 

These results are shown in figure 7 where it is seen that, for example, 
a 1-centimeter-thick bottom plate is worth about 0.3 centimeter of solu- 
tion height for this particular reactor. It is expected that the bottom 
plate worths for other values of R would be nearly the same as for the 
values for R = 475 since the fraction and spectrum of neutrons leaking 
from these do not differ greatly. It is estimated that the reported 
crit ical height for the NASA experiments could be too low by a maximum 
of about 0.5 centimeter since the wall thickness of the aluminum cylinder 
is 1/8 inch. Thus, because of the expected smallness of this correction, 
none of the Oak Ridge and NASA criticality data have been corrected for 
the criticality effect of the aluminum containment vessels. 


Since figure 6 is based on an energy- independent buckling, this 
implies an energy- independent extrapolation distance 5. This value of 
S can be deduced from experimental criticality data on unreflected 
cylinders if the top, bottom, and side extrapolation distances are 
assumed to be equal. Then, comparison with a wide range of reactor 
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dimensions and solution concentrations will establish the validity of a 
constant value for S. Figure 8 shows theoretical curves based on values 
of 5 ~ 2.5, 3.0, and 3.5 centimeters for unreflected 12-inch-diameter 
cylindrical reactors of UOgFg-H^O solution in which the critical height 

is plotted as a function of the parameter R. These curves are based on 
the following equation for the critical height h of the cylinders: 


where 



(35a) 


(35b) 


and r is the physical radius of the cylinders. For a given value of 
R, is evaluated from figure 6 for a K e pp of unity, and then, since 

the values of r and 5 have been fixed, the physical height h is 
given by equations (35a) and (35b). It is evident from figure 8 that 
these critical assemblies are quite sensitive to the value of the ex- 
trapolation distance which is used. For the 12 -inch-diameter cylindrical 
reactors considered the value of 5 = 3.0 centimeters seems to predict 
the experimental critical heights listed in table II quite well as a 
function of the atom ratio parameter R. This insensitivity of 5 with 
R is probably due to the fact that the hydrogen atom density, and hence 
the leakage fraction, changes but little over the range of R considered. 
Therefore, it was decided to use this value of 5 = 3.0 centimeters for 
the extrapolation distance for all values of the parameter R for the 
various diameters of the unreflected cylindrical critical assemblies of 
UO 2 F 2 -H 2 O. Reference 6 reports using values of 5 = 2.5 to 3.7 for the 
empirical kernel method of calculation. 


The curves of figure 9 represent calculations of the critical 
heights of cylindrical reactors for various diameters and as a function 
of the parameter R. These calculations are based on the theoretical 
buckling curve given as figure 6 and the assumed constant value of the 
extrapolation distance 6 = 3.0 centimeters. The experimental points, 
listed in table II, represent data from references 2, 3, 4, and the 
NASA solution reactor. The agreement for various diameter reactors as 
a function of R is quite good between the theoretical model described 
in this work and the various experiments. Thus, this provides the 
verification of the general buckling curve (fig. 6) and of the constant 
value of 3.0 centimeters for the extrapolation distance for these un- 
reflected cylindrical reactors. 
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Calculation of Nonthermal Fission Fractions, Leakage Fractions, 
Neutron Age, and Two- and Three-Group Macroscopic Cross 
Sections as a Function of Solution Concentration 
for UOgFg-H^O Assemblies 

In solving for the critical dimensions of unreflected UOgFg-E^O 

reactors using the consistent PI equations, the lethargy-dependent flux 
cp(u) and current j(u) are generated. Since all the information about 
a given reactor is contained by the flux and current, it is possible to 
derive many useful quantities pertaining to the UC^Fg-HgO reactors. 

Figures 10, 11, and 12 present some of this information. 

In figure 10 is plotted the nonthermal fission fraction as a func- 
tion of the parameter R for various cutoff energies. For a heavily 
loaded UOgFg-HgO system having R - 25, about 65 percent of the non- 
thermal fissions occur between energies of 10 7 ev to 0.125 ev, while 
about 43 percent of the total fissions occur between energies of 10 7 ev 
to 1.44 ev. For a UO 2 F 2 -H 2 O reactor having R = 700 the corresponding 

numbers are about 6 and 3 percent. 

Figure 11 is a plot of the fraction of the neutrons which escape 
from critical U02F2“R2° systems as a function of R. This fraction is 

nearly constant over the range of R from 25 to about 300, varying be- 
tween 44 and 47 percent. For larger values of R this function de- 
creases and is about 33 percent at R = 700. 

Figure 12 is a plot of the fission neutron age t as a function 
of the parameter R both for the indium resonance detector energy of 
1.44 ev and the thermal cutoff energy of 0.125 ev. From about R = 100 
to R = 700 the age is nearly constant. For values of R less than 
100 the age increases because of the decreased hydrogen atom concen- 
tration. In the limit as R -> co the 1. 44-ev curve would approach the 
fission neutron age in water. 

The neutron fluxes generated in solving for the critical dimensions 
of unreflected UO 2 F 2 --H 2 O systems can also be used as a weighting function 

for obtaining two- and three-group macroscopic cross-section parameters. 
These two- and three-group cross sections can be used in calculating 
reflected solution reactors. 


( 
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Both three- and two-group constants were generated as a function of 
R. Group I of the three-group split extends from u = 0 to u = 6.1, 
group II from u = 6.1 to u = 18.2, and the thermal group is centered 
about the lethargy of 19.8 corresponding to the room temperature kT of 
0.025 ev. For the two-group split, the thermal group is the same as 
that of the three-group split, and the fast group (group i) extends from 
lethargy u = 0 to u = 18.2. For the few-group calculations performed 
with a one -dimensional multigroup IBM 704 code, the parameters needed 
are the group diffusion coefficient iP, group macroscopic absorption 
cross section Egj, group macroscopic production cross section v5£, and 

the group removal cross section Z^* These group parameters are defined 
over the fluxes and currents generated for a given value of R as 



(36a) 


(36b) 


v (u)E^(u)cp(u)du 

( 36c ) 

Z cp(u)du 

w 0U P 



Q 

The quantity Z^ is evaluated by making a neutron balance over the 

group in question. Over the range of buckling considered, the group 
constants were not sensitive to the choice of buckling B^ 


The few-group constants obtained by the method described are shown 
in figures 13, 14, 15, and 16 for both the three- and two-group splits 
over the range of R values. The diffusion coefficient D is plotted 
in figure 13 and is relatively constant over the entire range of R for 
most of the fast groups; however, the effect of uranium-235 absorption 
is seen for all the R values for the thermal neutron group. The ab- 
sorption cross section is shown in figure 14 and varies strongly 

with the atom ratio R. The production cross section vZj is shown in 

figure 15 and exhibits similar behavior as the absorption cross section. 
The removal cross section Z^ is shown in figure 16 and increases as 
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the parameter R increases. Table III gives the corresponding three- 
and two-group data for water at room temperature. 


Comparison of Calculated and Experimental Critical Dimensions 
for Water-Reflected UOgPg-HgO Reactors 

In order to perform theoretical calculations on fully reflected 
cylinders, the reflector savings must he established as a function of 
reflector thickness and concentration parameter R. Using the three- 
group parameters plotted in figures 13 to 16 for various values of R 
and the water constants listed in table III, the critical core radius r 
was calculated for spherical reactors of solutions with various 

thicknesses of water using an IBM 704 diffusion code. The reflector 
savings RS is defined as 

RS - - r 

f BS 

2 

where for a given value of R the value of is obtained from fig- 

ure 6. The results of these calculations are shown in figure 17. The 
curves show a slight increase in the reflector savings with R. Por the 
completely reflected case (60 cm of water) the value of reflector savings 
is 6.5 centimeters at R - 25 and 6.0 centimeters at R = 700. Por the 
range of R values considered the reflector savings reached its maximum 
value at about 15 centimeters thickness of water. 



The critical heights h of fully reflected cylindrical assemblies 
were computed for cores of radius r as a function of the parameter R 
using the critical buckling curve of figure 6 and the infinite water 
reflector savings curve of figure 17. These results are shown in fig- 
ure 18 and were calculated with the following two formulas: 


B 


2 

z 



(38a) 

(38b) 


Experimentally obtained data listed in table 17 are also shown on the 
curves. The agreement is generally good. The dashed curves were com- 
puted by arbitrarily increasing the reflector savings by 0.5 centimeter 
for all values of R. The two sets of curves bracket the experimental 
data quite well and serve to indicate the sensitivity of the calculations 
with small changes in the value of the reflector savings. 
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CONCLUSION 

The lethargy-dependent equations of the consistent PI approxima- 
tion to the Boltzmann transport equation for slowing down neutrons, as 
developed in this present study, have been used as the basis of an IBM 
704 computer program. This program has been used to compute the age of 
fission neutrons in water as 26.5 square centimeters as compared to a 
recent experimental value given as 27.86 square centimeters. The 
Fourier transform of the slowing-down kernel for water to indium reso- 
nance energy calculated here compares well with the Fourier transform of 
the kernel for water measured by Hill, Roberts, and Fitch. This method 
of calculation has also been applied to uranyl fluoride - water solution 
critical assemblies. Theoretical results established for both unre- 
flected and fully reflected critical assemblies compare well with avail- 
able experimental data. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 22, 1961 
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APPENDIX - SYMBOLS 


A 


^(u) 


B 


2 

z 


D(u) 


D 


G 


dr 

du 

dv 


ratio of nuclear mass to neutron mass 

inelastic scattering parameter defined by eq. (33) 

lethargy-dependent material buckling defined by eqs. (8) 
and (9) 

transverse buckling for cylindrical geometry 

diffusion coefficient at lethargy u defined by eq. (22), 
cm 

neutron diffusion coefficient for group G defined by 
eq. (36a), cm 

spatial volume element 

element of lethargy 

velocity volume element 


E 


neutron energy 


E T 


P 


thermal cutoff energy separating thermal energy range 
from fast energy range 

spatial function defined by eq. (8b) 


f(u)du fraction of source neutrons in lethargy interval du 

about lethargy u 


f^(u) average cosine of angle between initial and final direc- 

tion of motion, measured in the center of mass coordi- 
nate system, for elastically scattered neutrons of 
lethargy u 

G spatial function defined by eq. (8a) 


H(v r v)drdv probability that a neutron in phase volume element 

(dr dv 1 ) will, after an inelastic scattering collision, 
be in the phase volume (dr dv) 

h physical height of cylindrical reactor 

j(u) lethargy-dependent neutron current defined by eq. (17b) 
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Fourier transform of neutron slowing-down kernel 
static criticality factor for a reactor 
Boltzmann constant 


K (%) 

Keff 

k 


N concentration of an element, atoms/cm 3 

n ( r , v ) dr dv probable number of neutrons in space volume element dr 

about position r and in velocity element dv about 
velocity v per unit time 


P(u) 


anisotropic slowing-down density defined by eq. (l9b) 


P(v’ ~>v)dr dv probability that a neutron in the phase volume element 

(dr dv’ ) will, after an elastic scattering collision, 
be in the phase volume (dr dv) 


P z (u-u r ) 


P Z (H) 

Q(r,u,H) 


lethargy-dependent coefficient of index 2, in expansion 
of elastic scattering kernel into a series of Legendre 
polynomials, defined by eq. (13) 

Legendre polynomial or argument p and index 2 

slowing-down density defined as number of elastic scat- 
tered neutrons which pass lethargy mark u per unit 
time per unit phase volume (dr dn) 


Q(u) isotropic slowing-down density defined by eq. (l9a) 

Q Z (Z, u) coefficient of index 2 in expansion of slowing-down 

density into a series of Legendre polynomials defined 
by eq. (7) 


R 

RS 

r 

— > 
r 

S(r, v)dr dv 


ratio of hydrogen atom concentration to uranium atom con- 
centration for the UOgPg-HgO solution reactors 

reflector savings defined by eq. (37) 

physical radius of cylindrical reactors 

position vector 

source of neutrons due to fission which are produced in 
phase volume element (dr dv) 

°K 


T 


moderator temperature. 
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contraction for lethargy difference u-u f 


neutron lethargy defined in terms of neutron energy 
by equation u = ln(E 0 /E), Eq usually taken as 10 7 


E 

ev 


cutoff lethargy corresponding to cutoff energy Erp 
neutron speed corresponding to energy E 


neutron velocity 


coordinate axis 


quantity related to mean squared lethargy increase per 
collision, defined by eq, (20b) 

Dirac delta function used in eq. (12) 

extrapolation distance defined by eqs. (35a) and (35b) 

Kronecker delta defined by eq. (5) 

reciprocal nuclear temperature defined by eq. (33) 

neutron transport mean free path defined by eq. (30) 


cosine of angle between Z-axis and direction of motion 
of neutron in planar geometry 

cosine of angle between initial and final directions of 
motion measured in laboratory coordinate system of a 
neutron undergoing an elastic scattering collision 

average cosine, measured in laboratory coordinate system, 
for elastically scattered neutrons 

mean of the product of lethargy increase in a collision 
times angle of scattering, defined by eq. (20c) 

average number of neutrons produced per fission occurring 
at energy E 

average increase in lethargy per collision 

quantity related to mean of the product of squared 
lethargy increase per collision times angle of scat- 
tering, defined by eq. (20d) 
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,G 


■'G 


macroscopic neutron absorption cross section, cm"^ 

macroscopic neutron elastic scattering cross section, cirT^ 

macroscopic neutron fission cross section, cm“^ 

macroscopic neutron inelastic scattering cross section, 
cm”l 

macroscopic neutron total cross section, E„ + E~ 0 + E-r^, 
cm"* 1 

macroscopic neutron absorption cross section for group G, 
defined by eq. (36b), cm“^ 

macroscopic neutron fission cross section for group G, 
defined by eq. (36c), cnf 1 

macroscopic neutron removal cross section for group G, 
cm"l 


a in 

t(E t , E) 


9 

y 

Q 


microscopic neutron absorption cross section, cm 2 

microscopic neutron elastic scattering cross section, cm 2 

microscopic neutron fission cross section, cm 2 

microscopic neutron inelastic scattering process, cm 2 

neutron age, defined as one-sixth the mean squared 
slowing -down distance to energy E of a neutron 
originating from a point source emitting neutrons of 
energy E f in an infinite medium 

neutron flux, nv 

unit vector in direction of motion of neutron 
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TABLE I. - FOURIER TRANSFORM K(B m ) OF SLOWING -DOWN 
KERNEL FOR FISSION NEUTRONS IN WATER AND IN 
U°2F2~H 2 0 SOLUTIONS WITH ATOM RATIOS 
OF HYDROGEN TO URANIUM-235, R, 


EQUAL TO 50, 150 



K(Bj 

H 2 0(R=oo) 

R = 150 

R = 50 

0 

1.0 

1.0 

1.0 

.03142 

.9745 

.9748 

.9752 

.06283 

.9047 

.9057 

.9069 

.12566 

.6944 

.6970 

.7000 

.20944 

.4188 

.4220 

. 4263 

.26180 

.2881 

.2922 

.2970 

. 31416 

.1923 

.1963 

.2013 

.3927 

.0995 

.1030 

.1079 

. 5236 

.0285 

.0310 

.0349 

.7854 

.0005 

.0013 

.0028 




TABLE II. - EXPERIMEHTAL CRITICALITY DATA FOR U0 2 F 2 ~H 2 0 URREFLECTED 


SOLOTIOF REACTORS CORTAUJED BY ALOMINIM VESSEIS 


Cylinder 

diameter, 

in. 

Atom ratio, 

R - R 2 /^ 235 

Critical 

height, 

cm 

Reference 

Cylinder 

diameter, 

in. 

Atom ratio, 

R - FV - 235 

Critical 

height, 

cm 

Reference 

10 

52.9 

34.0 

2, p. 75 

15 

169.0 

18.5 

2, p. 75 


169.0 

41.2 

2, p. 75 


328.7 

21.7 

1 







499.0 

27.4 

i 


27.1 

38.9 

3, p. 24 


755.0 

43.6 

7 


44.3 

35.1 


1 






73.4 

33.7 


f 


27.1 

18.5 

3, P. 25 







44.3 

17.9 

1 


43.2 

a 34,9 

4, p. 

73-74 


74.6 

16.8 

Y 


44.7 

a 34.7 








50.1 

34.8 




50.1 

17.9 

4, p. 73 


51.5 

a 33. 5 




331 

22.9 

4, p. 73 


55.4 

34.3 








60.8 

34.1 



20 

27.1 

15.8 

3, p. 25 


66.1 

34.1 




44.3 

15.0 

1 


71.5 

34.1 




73.4 

15.2 

Y 


83.1 

34.4 








85.7 

34.9 




50.1 

15.4 

4, p. 73 







60.8 

15.3 

1 

12 

44.3 

23.2 

3, I 

. 25 


325 

18.7 



50.1 

22.6 

4, I 

. 73 

30 

44.3 

13.7 

3, p. 25 


51.5 

22.6 




72.4 

13.9 

3, p. 25 


55.4 

a 22. 6 








60.8 

22.7 




50.1 

13.8 

4, p. 73 


331 

32.8 

* 


51,5 

a 13. 3 

[ 







331 

16.3 

Y 


185 

25.3 

Present 






218 

26.6 

work 






250 

28.3 








253 

28.5 








285 

30.2 








320 

32.0 








365 

35.4 








405 

37.0 








435 

43.6 








455 

46.2 








460 

47.0 








465 

48.2 








485 

48.8 








^Value too low according to ref. 4. 



TABLE III* - THREE- AND TWO -GROUP CONSTANTS FOR WATER REFLECTOR 


Lethargy, 

u 

Diffusion 

coefficient, 

D P 

Ah sorption 
cross section, 

% 

Removal 
cross section, 
vG 

0-18.2 

1.274 

0.00066 

0.04598 

0-6.1 

1.722 

.00012 

. 07837 

6.1-18.2 

.612 

.00145 

. 11392 

19.8 

.1329 

.01962 



TABLE IV. - EXPERIMENTAL CRITICALITY DATA FOR U0 2 F 2 -H 2 0 REACTORS 
FULLY REFLECTED BY AN INFINITE WATER REFLECTOR 


Cylinder 

diameter, 

in. 

Atom ratio, 

R « N^rf 1-235 

Critical 

height, 

cm 

Reference 

8 

26.2 

21.5 

2, p, 69 


29.9 

20.7 




52.9 

19.5 




58.8 

20.5 




99.5 

22.4 




192.0 

2Q.1 

t 



290.0 

40.1 




43.2 

18.6 

4, p. 72 

10 

52.9 

13.4 

2. p. 69 


328.7 

22. 4 

J 



499 

35.2 

\ 



43.2 

12.5 

4, I 

). 72 


51.5 

11.4 

4, p. 72 

15 

52.9 

7.90 

2, p 

). 69 


56.7 

8.50 




221.0 

11.30 




499.0 

16.90 




755.0 

27.10 

T 
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Figure 3. - Uranium-235 atom concentration in U0 2 F2 -H 2° solutions. Uranium enriched to 0.932 uranium-235; water at 25° C; density, 0.9970 gram 
per cubic centimeter. 
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Figure 5. - Lethargy-dependent flux <p(u) plotted against lethargy u for a critical U0 2 F 2 solution reactor of N H/ N U-235 „ 472 . Source normalization, 1 fission 

/•U 

neutron. Thermal neutron flux ©m « 5.661; / cp(u)du = 14.333. 
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Figure 8. - Critical height for several values of extrapolation length 

plotted against atom ratio N H /N U ~ 235 for 12-inch-diameter unreflected 
U0 2 F2 solution reactors. Theoretical curves based on total buckling 

from figure 6 . 



40 



Figure 9. - Critical height plotted against atom ratio N h /N u ~ 235 for various 
diameters of unreflected UO 2 F 2 solution reactors. (Theoretical curves based 
on total buckling from figure 6 and extrapolation distance 6 = 3.0 cm.) 









Nonthermal fission fraction . 


E-690 



Figure 10. - Fraction of nonthermal fissions to total fissions plotted against atom ratio 
FH/nU- 235 for several cutoff energies for U02F2~H20 solution reactors. 
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Figure 12. - Fission neutron age t as a function of atom ratio to indium resonance 

energy (1.44 ev) and to thermal cutoff energy (0.125 ev = 5 kT) for U0?F?-H?0 solution 
reactors. 
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Absorption cross section. 
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Figure 14. - Two and three-group absorption cross sections 2 a plotted against atom ratio 
for UO 2 F 2 solution reactor. 
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Figure 15. - Two- and three-group production cross sections v2f> plotted 
against atom ratio N h /N u - 235 for U0 2 F 2 solution reactors 


t 


< 


% 



• uocpBS .1 uoxq.nios 2^on *ioj g£ 2 - n H/ H M 

OTq.'ej uiogB jo uoxq.oury: a sb 2 suoxpoas ssouo um.op-Sutmo.xs dnouS-asuqg puu -o.mx - -gp axnSxjj 


S22-Il M /l# ‘°?V* 











Critical height, cm 


O Ref. 2 

A Ref. 4 

Computed using theoretical 

reflector savings data 
of fig. 17 

Computed by increasing 

reflector savings data 
by 0.5 cm 


Diameter, 

in. 


O 

-A 




C 40 60 80 100 200 400 600 8 

Atom ratio, N^/n^ -233 

Figure 18. - Critical height plotted against atom ratio N H /N U-235 for various 
diameters of cylindrical UO 2 F 2 solution reactors having both end and side 
reflectors of water (infinite thickness). Theoretical curves obtained using 
calculated values of total buckling given in figure 6 and reflector savings 
given in figure 17. 
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